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ABSTRACT 

Near-field acoustical holography (NAH) requires the measurement of the 
near-field pressure field over a conformal and closed surface in order to 
recover the acoustic field on a nearby surface. We are interested in the 
reconstruction of the acoustic field over the fuselage of a Boeing 757 
airplane when pressure data is available over an array of microphones 
that are conformal to the fuselage surface. In this case the strict NAH 
theory does not hold, but still there are techniques used to overcome 
this difficulty. The best known is patch NAH, which has been used for 
planar surfaces. In this work we will discuss two new techniques used 
for surfaces with an arbitrarily shape: patch inverse boundary element 
methods (IBEM) and patch equivalent sources method (ESM). We will 
discuss the theoretical justification of the method and show reconstruc- 
tions for in-flight data taken inside a Boeing 757 airplane. 

1 INTRODUCTION 

Denote as T the surface of a cylinder as shown in figure 1. For a time-harmonic (e~ luJt ) dis- 
turbance of frequency u the sound pressure p satisfies the homogeneous Helmholtz equation 
inside the fuselage 

A p + k 2 p = 0, (1) 

where k = u/c is the wave number and c the constant for the speed of sound. 

The conventional solution of interior NAH will require that the pressure should be mea- 
sured on a closed surface interior to F. Ideally this measurement surface will be very close to 
T in order to obtain a high-resolution reconstruction. However, practical implementations 
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Figure 1: Cylindrical surface T and the patch surface T p . 


does not allow a closed and all encompassing measurement surface. Thus a smaller mea- 
surement surface, or “patch” surface, is used instead. In this work we will consider the case 
when measurements are taken over a patch To, where this patch is conformal the surface F p 
(see figure 1). 

We use the integral representation [1, 2, 3] for a point x = (xi,X 2 ,Xs) E Tq 


P(x) = / d>(x,yV(y)dS(y), 
Jtp 


where 

<h(x,y) = 

Similarly the normal velocity is given by 

5$(x,y) 


exp {ik |x — y| 
47T |x — y| 


ipujv^x) = 

Jr 


rp aa(x) ^ (y) ^ (y) + 2^ (X) ’ XerP ’ 


( 2 ) 


( 3 ) 


where n is the outward unit normal and p the constant of the medium density. 

A solution p of (1) can be represented by the method of equivalent sources[4, 5]. In this 
method we are required to define a source surface T s that lies at a constant distance S from 
T p as shown in figure 2. This method uses the representation 

N s 

K x ) = Z J)’ z i e ( 4 ) 

3 = 1 


where N s > 0 is the number of sources and qj,j = 1 ,...,N S the source coefficients. The 
normal derivative is calculated using 


N s 

ipuv(y) = Qj 
3=1 
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2 NUMERICAL DISCRETIZATION 

The boundary surface F p is decomposed into quadrilateral elements with four nodes. The 
boundary element method with iso-parametric linear functions[3] is selected for interpolating 
the geometric and acoustical quantities. Given M pressure measurements on T 0 , represented 
as p, recover N pressure and normal velocity points on F p . represented as p s and v s respec- 
tively. When x G Tq, (2) gives the matrix equation 

[S] P> = P (6) 

where [S] is a M x N complex matrix and ip is the column vector of N entries that represent 
values of the density <p on F. Similarly, when x G F, (3) produce the matrix equation 

V* = [K] v , (7) 

where [K] is an N x N complex matrix. 

The source surface T s does not need to be decomposed into surface elements, instead we 
just need to create a grid with N s points distributed over T s . Then is simple to obtain the 
following matrix system 

[G] q = p, ( 8 ) 

where the coefficients of the M x N s complex matrix [G] are given by 

Gij = $(xi,Zi). 

Here are the points in T 0 and z j are the points in T s . Similarly we obtain the relationship 

v s = — [G s?; ] q (9) 
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where the coefficients of the N x N s complex matrix [G 6 ""] are given by 


r sv = 

13 dn(yi) 


and y, are the points in T 

3 NUMERICAL REGULARIZATION 

For the experimental problem, the exact pressure p is perturbed by measurement errors. We 
denote the measured pressure as p m . If the elements of the perturbation e = p m — p are 
Gaussian (unbiased and uncorrelated) with covariance matrix a q [I], then E ( 1 1 e 1 1 1 
(here 1 1 . 1 1 2 is the 2-norm). It is well known that the linear systems on (6) and (8) are ill-posed, 
i.e., the errors in p m will be amplified in the solution <p or q, and in most of the cases the 
recovery will be useless. 

To avoid the amplification of the measurement errors, special regularization methods are 
used to find the solution of these linear systems. The best known implementation of these 
methods requires the use of the singular value decomposition (SVD) 

[S] = [U 1 ][S 1 ][V 1 ] H , [g] = [u 2 ] [e 2 ] [v 2 ] H 

where [U 1,2 ], [V 1,2 ] are unitary matrices and [£ 1,2 ] is a diagonal matrix containing the singu- 
lar values G{ in order of non-decresing magnitude. Regularization methods are implemented 
using the explicit formula 

^ = [V 1 ] [F 1,Q ] [S 1 ] [U 1 ]" , q a = [V 2 ] [F 2 ’"] [S 2 ] [U 2 ]" , (10) 

Here the diagonal matrices [F 1 ,a ],[F 2,a ] contain the filter factors[6] which are used to reduce 
the effect of the measurement errors in the reconstruction. The parameter a > 0 is called the 
regularization parameter, and should be chosen correctly. There are several types of filter 
factors used for different ill-posed problems, but in this work we will use the filter factors of 
Tikhonov with a high-pass filter [7]. 

4 PHYSICAL EXPERIMENTS 
4.1 Tin can 

The experimental configuration for the holographic measurement is similar to the previous 
work of Herdic[8]. The surface T is an aluminium stiffened cylindrical shell (0.81m radius, 
2.55m length and the shell thickness varies between 0.8 and 1.2mm) excited by a point 
force applied to a rib/stringer intersection at one end of the cylinder. The measurements 
were conducted using a chirp waveform over a band from 10 to 1000Hz with 0.61 Hz resolu- 
tion. For this interior NAH problem, the measurement surface To is a cylindrical array of 
0.7045m radius, 123.75 degrees angle and 0.74m length as shown in figure 2. The pressure 
measurements on r 0 are a grid of 12 points over the radius and 10 points over the length. 
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We use this experiment for the validation of the patch IBEM and ESM methods, since 
the normal velocity is known from vibrometer measurements at the surface of the cylinder. 

is a patch of the surface F that is directly in front of r 0 . F p is a cylindrical surface of 
0.81m radius, with the same angle and length than r 0 . The normal velocity is measured 
over 120 points in distributed as in r 0 . We will consider the use of extension points N x 
in T p (and F® for ESM), which will make the cylindrical surface angle 146.25, 168.75 degrees 
and respectively the length 0.91, 1.08m for N x = 1,2 (see figure 2). 

The vector v s contain the measured normal velocity on T p with N x = 0. Notice that 
the previous methods will provide the reconstructed normal velocity v® for (even over 
extended points N x > 0), but the relative error 1 1 v® — v® 1 1 2 / 1 1 v® 1 1 2 x 100 will be considered over 
the part of that was not extended (N x = 0). Since we know the exact v®, we choose in all 
our reconstructions the optimal regularization parameter a for the Tikhonov regularization 
with a high-pass filter in (10). 



Figure 3: Patch IBEM reconstruction for different frequencies using extension points N x . 


In figure 3 we show the relative error of reconstructed normal velocity using patch IBEM 
for 44 frequencies equally spaced from 88.5Hz to 180Hz. This figure clearly shows the 
importance of the extension points. When N x = 0 the errors are bigger for all frequencies. 
N x = 1,2 have similar errors, but N x = 2 does slightly better. If we utilize more extension 
points we will observe that the errors will not decrease, but instead the columns of the matrix 
system in (7) will increase. The increase in the columns of the matrix [S] in (7) will increase 
considerably the memory requirements of patch IBEM. For that reason is just necessary to 
use N x = 2 with patch IBEM. 

In figure 4 we show the relative error of the reconstructed normal velocity using patch 
ESM with 8 = 0.08m for 44 frequencies equally spaced from 88.5Hz to 180Hz. This figure 
shows that the reconstruction errors are similar to the reconstruction errors of patch IBEM 
for different extension points. When N x = 0 the errors are higher for all frequencies, and 
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Figure 4'- Patch ESM with 8 = 0.08m reconstruction for different frequencies using extension 
points N x . 


for N x > 0 the errors reduce. Again as for the patch IBEM method, there will not be a 
considerable decrease in the reconstruction point for N x greater than 2. 

In figure 5 we show the comparison between the errors of patch IBEM and ESM with 
8 = 0.08m. In both methods we use N x = 2 and we found that they produce similar errors 
for all frequencies. 

4.2 Boeing 757 fuselage 

The surface T is the fuselage (1.74m radius) excited by a acoustic sound source. The mea- 
surements are available at a cylindrical array of 1.73m radius, 23 degrees angle and 0.27m 
length (as the cylindrical array in figure 2). The pressure measurements on T 0 are a grid of 
12 points over the radius and 10 points over the length. 

T p is a patch of the surface T that is directly in front of T 0 . T p is a cylindrical surface 
of 1.74m radius, with the same angle and length than To- As in the previous experiment we 
consider the use of extension points N x in T p (and P for ESM). 

In figure 6 we show the normal velocity reconstruction v;l in P’ using patch IBEM with 
extension points N x = 2. Finally in figure 7 we show the normal velocity reconstruction v;l 
in T p using ESM with extension points N x = 2 and 8 = 0.0254m. 


5 CONCLUSION 

In this work we have shown two methods that are used for the reconstruction of the normal 
velocity when measurements are available over a patch: patch IBEM and ESM. It was found 
that both methods produce similar reconstruction errors, and also in both methods these 
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Figure 5: patch IB EM and ESM reconstruction for different frequencies using extension point 
N x = 2. 


errors can be reduced by the use of extension points in the reconstruction surface P’ (and 
source surface P for ESM). 

Patch ESM does not require the time-expensive calculations of the integral as in IBEM. 
For that reason we have shown that patch ESM can be numerically more attractive than 
patch IBEM, when patch measurements are available. 
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Figure 6: Normal Velocity reconstruction using patch IBEM with extension points N x = 2 
for some selected frequencies. 
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Figure 7: Normal Velocity reconstruction using patch ESM with extension points N x = 2 and 
8 = 0.0254m for some selected frequencies. 
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